IMAGE PROCESSING METHOD , IMAGE PROCESSOR, AND STORAGE MEDIUM 

THEREOF 

BACKGROUND OF THE INVENTION 

1 . Field of the Invention 

The present invention relates to an image processing method, 
an image processor and a storage medium to perform half tone 
processing on input images using an error diffusion method, and 
more particularly to an image processing method, an image 
processor and a storage medium thereof wherein the geometric 
texture and graininess of output images are improved with the 
error diffusion method. 

2 . Description of the Related Art 

Recently full color image processing has become possible 
with advances in color image processing technology. For example, 
processing with 2 56 grayscale levels is performed for each one 
of cyan C, magenta M and yellow Y. However, such a device as 
a display and a printer, to display and print these colors, cannot 
display or print full color images since the number of grayscale 
levels is limited. Therefore, such a device as a display and 
a printer, where the number of grayscale levels which can be 
reproduced is limited, converts a full color input image to a 
half tone image by the error diffusion method. 



1 




The error diffusion method uses an algorithm to obtain half 
tone images where grayscale is satisfactorily preserved by 
diffusing the quantization errors , which are generated when input 
images are quantized by the number of grayscale levels possible 
5 in the output device, into the peripheral pixels. 

However, a shortcoming of the error diffusion method is 
the generation of geometric texture , such as the texture generated 
in the low density dot area (area where the density of block 
dots is low) of the output image (this texture is called "worm" ) , 

10 and the texture generated in the high density dot area (area 
where the density of black dots is high) (this texture is called 
"fingerprint"), and many improvements have been attempted in 
the prior art. The improvement techniques which have been 
proposed will be described below. 

15 The first method is a two-way scan (Ulichney, 1987) which 

scans an input image in two ways, where worm is decreased since 
the diffusion direction of the quantization errors is dispersed 
in the low dot density area. 

The second method is an error filter switching system, which 

20 switches the error filters to diffuse the quantization errors 
into the peripheral pixels (R. Eschbach, "Reduction of artifacts 
in error diffusion by means of input dependent areas", Journal 
of Electronic Imaging, 2 (4), pp. 352 - 358, (1993)). Depending 
on the input pixels and the peripheral pixels thereof, several 

25 types of error filters are used to decrease worm and fingerprint. 
Worm decreases as the filter to be used becomes larger in size. 
A system combining the two-way scan system and the error 
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filter switching system has also been proposed (e.g. Japanese 
Patent Laid-Open Application No. 9-65126). 

The third method is a sine wave superimposing system (Y. 
Kishimoto, M. Nose, R. Saito and H. Kotera, "Improved error 
5 diffusion method with AM/FM periodic noise", Rroc. IS&T SNIP15, 
pp. 366 - 369 (1999)), where dots are corrected and arranged 
by superimposing a sine wave for modulating the amplitude and 
the frequency on the input pixel values before quantization 
processing. The volume of calculation is relatively low. 
10 The above mentioned prior art, however, has the following 

problems - 

In the first method, two-way scanning, worm is decreased 
since the diffusion direction of quantization errors is dispersed 
in the low dot density area, but fingerprint is amplified by 
15 the interference of quantization errors in the high dot density 
area . 

In the second method, which is a combination of the error 

filter switching system and the two-way scan, geometric texture 

(worm and fingerprint ) decreases but uniform cyclic dot patterns , 
20 such as half tone dots, cannot be created in the low dot density 

area, the graininess of dots becomes salient, and the color 

reproduction range does not improve. 

In the third method, the sine wave superimposing system, 

sine waves are applied simply according to the intensity of the 
25 input level, so the edges of the output image become unclear, 

and images where density varies greatly cannot be improved very 

much. 
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With the above mentioned techniques alone, it is difficult 
to improve both worm and fingerprint simultaneously, and the 
indiscriminate combination of techniques merely increases the 
volume of calculation. 

5 

SUMMARY OF THE INVENTION 



It is an object of the present invention to provide an image 
processing method, an image processor and a storage medium not 
10 only to prevent geometric texture, but also to prevent graininess 
from deteriorating. 

It is another object of the present invention to provide 
an image processing method, an image processor and a storage 
medium to obtain images with low graininess , close to half tone 
15 dots, in a low dot density area. 

It is still another ob j ect of the present invention to provide 
an image processing method, an image processor, and a storage 
medium to prevent unclear edges of an output image and to present 
the effect of error diffusion even for images with highly varied 
20 density. 

To achieve these objects, an image processing method to 
perform half tone processing on input images using an error 
diffusion method according to the present invention comprises 
steps of: generating a sine wave where the amplitude and the 
25 frequency are modulated according to the average value of the 
target pixel value and the peripheral pixel values; adding the 
diffused quantization errors, the target pixel value and the 
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sine wave; quantizing the addition result by a predetermined 
number of grayscale levels; and calculating the quantization 
errors of the peripheral pixels from the errors by the 
quantization . 

5 An image processor which performs half tone processing on 

input images using an error diffusion method according to the 
present invention comprises a memory for storing the input image 
and a processing part for performing the half tone processing 
on the input images , wherein the processing part adds a sine 

10 wave where the amplitude and the frequency are modulated according 
to the average value of the target pixel value and the peripheral 
pixel values, the diffused quantization errors, and the target 
pixel value, then quantizes the addition result by a predetermined 
number of grayscale levels, and calculates the quantization 

15 errors of the peripheral pixels from the errors by the 
quantization. 

A program of a storage medium for storing a program to perform 
half tone processing on input images using an error diffusion 
method according to the present invention comprises: a program 

20 for generating a sine wave where the amplitude and the frequency 
are modulated according to the average value of the target pixel 
value and the peripheral pixel values; a program for adding the 
diffused quantization errors, the target pixel value and the 
sine wave; a program for quantizing the addition result by a 

25 predetermined number of grayscale levels; and a program for 
calculating the quantization errors of the peripheral pixels 
from the errors by the quantization. 
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In the present invention , the amplitude and the frequency 
of a sine wave to be superimposed are modulated using an average 
value of the target pixel value and the peripheral pixel values. 
Since the average value indicates the level of dot density, a 
5 sine wave with amplitude and frequency according to the dot density 
can be superimposed. Therefore, where error diffusion is 
performed to maintain the grayscale levels, uniform and cyclic 
dot patterns, as in half tone dots, can be arranged in the low 
dot density area, and a high frequency dot pattern, similar to 
10 a checkered pattern, can be arranged in the high dot density 
area, and as a result, a visually better half tone image can 
be obtained. 

In the image processing method according to the present 
invention, the step of generating a sine wave may further comprise 

15 a step of calculating the average value of the target pixel value 
and the peripheral pixel values, a step of calculating the 
difference between the target pixel value and the average value, 
and a step of modulating the amplitude of the sine wave by the 
average value and the difference. Since difference is used, the 

20 edges of the image can be detected and the effect of the sine 
wave can be decreased so that a forced arrangement of dots at 
the edges can be prevented, and clear edges can be stored. As 
a result, a good half tone image can be obtained, even for images 
with highly varied densities. 

25 In the image processing method according to the present 

invention, the step of generating a sine wave may further comprise 
a step of independently controlling the cycle in the horizontal 
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direction and the cycle in the vertical direction of the sine 
wave. Since the horizontal direction and the vertical direction 
are independently controlled , the sine wave functions in two 
dimensions , and the arrangement of dots can be controlled in 
5 two dimensions using the sine wave. As a result, a better half 
tone image can be obtained in a two-dimensional plane. 

In the image processing method according to the present 
invention, the step of generating a sine wave may further comprise 
a step of setting different initial phases of the sine wave for 

10 the input image of each color . When color images are manipulated, 
dots can be forcibly distributed so that the dots of each color 
/ do not overlap, and the dots of each color can be arranged at 
relatively equal intervals. As a result, the spatial frequency 
of a half tone image can be increased, graininess can be decreased, 

15 and the color gamut expanded. 

The image processing method of the present invention may 
further comprise a step of reading the target pixel value and 
the peripheral pixel values by scanning the input image in two 
ways, and a step of calculating the quantization errors of the 

20 peripheral pixels may further comprise a step of calculating 
the quantization errors using an error filter, which is selected 
from a plurality of error filters according to the target image 
value. In this aspect, the increase in calculation volume can 
be controlled and geometric texture can be controlled by using 

25 both two-way scanning and a small size error filter switching 
system. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

Fig. 1 is a block diagram depicting an image processor of 
an embodiment of the present invention; 
5 Fig. 2 is a functional block diagram depicting the image 

processing unit in Fig. 1; 

Fig. 3 is a diagram depicting the two-way scan in Fig. 2; 

Fig. 4 is a block diagram depicting the most important part 
in Fig. 2; 

O 10 Fig. 5A, 5B and 5C are diagrams depicting the dot density 

. -i 

; 4 of the present invention; 

p| Fig. 6A and 6B are diagrams depicting the relationship 

^ between the average value and the dot density; 

Fig. 7A and 7B are diagrams depicting the amplitude 
15 calculation of the AM— FM screen function; 
H Fig. 8A and 8B are diagrams depicting the operation of the 

3 amplitude calculation of the AM-FM screen function; 

Fig. 9A and 9B are diagrams depicting the amplitude control 
of the AM-FM screen function; 
20 Fig . 1 0 is a diagram depicting the operation of the amplitude 

control of the AM-FM screen function; 

Fig. 11 is a diagram depicting the frequency calculation 
of the AM-FM screen function; 

Fig. 12 is a diagram depicting the frequency calculation 
25 of the AM-FM screen function; 

Fig. 13 is a diagram depicting the two-dimensional control 
of the AM-FM screen function; 
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Fig. 14 is a diagram depicting the amplitude and the frequency 
control of the AM-FM screen function; 

Fig. 15 is an image processing flowchart of an embodiment 
of the present invention; 
5 Fig. 16 is a diagram depicting color image processing of 

an embodiment of the present invention; 

Fig. 17 is a diagram depicting the color image processing 
in Fig. 16; 

Fig. 18 is a diagram depicting an output pattern of a 
10 comparison example; 

Fig. 19 is a diagram depicting an output pattern of another 
comparison example; and 

Fig. 20 is a diagram depicting an output pattern of an 
embodiment of the present invention. 
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DESCRIPTION OF THE PREFERRED EMBODIMENTS 



The embodiments of the present invention will now be 
described in the sequence of image processing and the AM-FM screen 
20 function with reference to the accompanying drawings. 
[Image processing] 

Fig. 1 is a block diagram depicting an embodiment of the 
present invention, Fig. 2 is a functional block diagram of the 
image processor thereof, Fig. 3 is a diagram depicting two-way 
25 scanning, and Fig. 4 is a block diagram depicting a part of the 
image processor. 

Fig. 1 is an example when the present invention is applied 
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to a color ink jet printer, where the memory 20 stores the input 
image as a full color image having 0 - 255 density values (0: 
black, 255: white) in each C (cyan), M (magenta) and Y (yellow) 
plane. The image processing unit 21, which comprises an MPU, 
5 performs Q value quant inze processing, and the quantization 
threshold is set at a constant interval. For example, in the 
case of binary (base-2) processing, the threshold is set at the 
12 7.5 level, and in the case of sextic (base-6) processing, the 
threshold is set at the 25.5, 7 6.5, 127.5, 178.5 and 229.5 levels. 

10 The output memory 22 stores the half tone processed output 

image to the C, M and Y planes. The color ink jet head 23 jets 
ink drops in each color according to the output image, and prints 
to the printing medium, which is not illustrated. The MPU 21 
processes each color of the input image for each color. The 

15 memories 20 and 22 may be constituted as one memory. 

The image processing unit 21 will be described with reference 
to Fig. 2. In Fig. 2, '1' is an input image. '2' is an addition 
circuit for adding the input pixel value I and the quantization 
errors E determined through the error filter 6, '3' is a 

20 quantization circuit for quantizing the value after addition 
by the addition circuit 2 using quantization thresholds, '4' 
is a subtraction circuit for calculating the quantization errors , 
'5' is an error buffer for recording the quantization errors, 
' 6 ' is an error filter where an error diffusion ratio is recorded, 

25 ' 7 ' is a density level detection circuit for detecting the values 
of the target pixel and the peripheral pixels in the input image, 
'8' is a circuit for generating the AM-FM screen function (sine 
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wave) based on the value from the detection circuit 7 , and '9' 
is a buffer for storing the cycle of the sine wave. 

Before a description of the operation of this configuration , 
the two-way scanning used for this embodiment will be described. 
5 Two-way scanning is a scanning method of reciprocating scanning 
of each scanning line. The cause of worm generated in a low dot 
density area is mainly a deviation of the quantization error 
diffusion direction in a certain direction, but the use of this 
scanning method decreases worm since the dispersion of the error 

10 diffusion direction improves. 

When two-way scan is used, however, fingerprint may be 
amplified as a side effect, so two types of small size error 
filters, 60 and 61, shown in Fig. 4, are provided in the error 
filter 6. Fingerprint can be controlled by selecting and using 

15 the small size error filters 60 and 61 according to the value 
detected by the detection circuit 7 in Fig. 2, that is, according 
to the density level of the target image value. 

The error filter 60 in Fig. 4 has a characteristic which 
increases the uniformity of dots in a highlight area (low dot 

20 density area), and the error filter 61 has a characteristic to 
bring out the high frequency characteristic of dots in the mid-tone 
area (high dot density area) . As Fig. 4 shows, the low dot density 
area is an area where the number of black dots is low (number 
of white dots is high), and an area where the number of black 

25 dots is high (number of white dots is low), for which the error 
filter 60 is selected. The high dot density area is the area 
in between, where image changes considerably. The error filter 
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61 is selected for this area. 

By using both of the above two methods, two-way scan and 
small size error filter switching operation/ an increase in the 
calculation volume can be minimized and such problems as geometric 
5 texture as worm and fingerprint can be solved. 

Then during quantization processing, the AM-FM screen 
function (sine wave) having an amplitude and a frequency, which 
depends on the value obtained from the density level detection 
circuit 7 in Fig. 2, is calculated by the generation circuit 

10 8, and the calculated value along with the quantization errors 
are superimposed on the addition circuit 2. In this way, dot 
patterns with low graininess, similar to half tone dots, can 
be obtained in the low dot density area of an output image. 
Now the operation of the configuration in Fig. 2 will be 

15 described. At first, the scanning direction is determined. The 
first scanning direction can be in the right direction, as in 
a normal raster scan. Then the density level detection circuit 
7 detects the density levels of the target pixel value and the 
peripheral values of the input image. The error filter 60 or 

20 61 (6 in Fig. 2) to be used for the pixel value is determined 
according to the density level of the target pixel value detected 
by the circuit 7,. 

According to the error filter 60 or 61, the quantization 
errors stored by the error buffer 5 are weighted and added by 

25 the circuit 6. Here, the average density level Avg of the target 
pixel value detected by the detection circuit 7 and of the 
peripheral pixel values 8 thereof, is calculated by the low pass 
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filter 70 of the detection circuit 7 in Fig, 4 , and the difference 

6 between the target pixel value I and the average density level 
Avg is calculated by the difference circuit 71. 

Based on the average density level Avg and the difference 

5 level 6 obtained there, the generation circuit 8 calculates the 
amplitude and the frequency of the AM-FM screen function, and 
determines the density level to be superimposed. As mentioned 
later in Fig. 4, the amplitude calculation circuit 80 calculates 
the amplitude from the average density level Avg and the difference 

10 level 6, and the frequency calculation circuit 81 calculates the 
frequency from the average density level Avg. And the function 
calculation circuit 82 calculates the function to be superimposed 
from the amplitude and the frequency. 

In Fig. 2, the target pixel value I (x, y) , the quantization 

15 errors obtained from the circuit 6 , and the AM-FM screen function 
value calculated by the circuit 8, are added by the addition 
circuit 2 . The value added by the addition circuit 2 is quantized 
to the level determined by the threshold of the quantization 
circuit 3, and the value becomes the pixel value of the output 

20 image . 

The circuit 4 subtracts density level determined by the 
circuit 3 from the value determined by the circuit 2. The value 
calculated by the subtraction circuit 4 is recorded in the error 
buffer 5 . After this operation, the target pixel value is shifted 
25 to the right by 1 , and the above series of operations is repeated. 

When the above series of operations is repeated until the 
last pixel of the main scanning line, the scanning direction 
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of the scanning line to be processed next is reversed (in other 
words, from right to left) , and the same operation is performed. 
Since operation is in the opposite direction this time, the shapes 
of the error filters 60 and 61 in 6 are also reversed. 
5 [AM-FM screen function] 

The amplitude characteristic and the frequency 
characteristic of the AM-FM screen function will now be described. 
At first, the amplitude (hereafter Amp) characteristic of the 
AM-FM screen function will be described with reference to Fig. 

10 5 to Fig. 10. For the amplitude Amp, the average value (hereafter 
Avg) of the peripheral 9 pixels using the low pass filter 70 
in the detection circuit 7, and the amplitude Amp are controlled 
using the average value. The AM-FM screen function is set such 
that the amplitude is high in the low dot density area, and is 

15 low in the high dot density area. 

As Fig. 6A shows, binary processing where Q is "2" is used 
as an example for the description, the number of black dots is 
less as the average value becomes greater than the threshold, 
and this is defined as a low dot density area. Fig. 5A shows 

20 this state in a 3 x 3 matrix. The number of white dots is less 
as the average value becomes smaller than the threshold, and 
this is also defined as a low dot density area. This is shown 
in Fig. 5C. When the average value is near the threshold, black 
and white dots coexist, as shown in Fig. 5B. This is defined 

25 as a high dot density area. Fig. 6B shows the case of quaternary 
(base-4 ) processing. 

The amplitude of the AM-FM screen function indicates the 
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intensity to with which dots are corrected and arranged dots 
in a highlight part (low dot density area). If the amplitude 
of the sine wave is simply determined according to the intensity 
of the input level, the edges of the output image may become 

5 unclear. So the difference 6 between the target pixel value and 
the value calculated by the low pass filter 70 is used, and the 
difference is also set as the amplitude parameter of the sine 
wave. By controlling the amplitude in this way, clear edges are 
maintained and the dispersion of dots becomes more stable. 
10 The amplitude Amp is given by the following formula. 

Amp = C(constant) • (mod(Avg, 255/ (Q-l ) )-255/2 (Q-l ) ) y 
/(6 4 +l) (1) 

Here mod (a, b) is a remainder when 'a' is divided by *b'. 
This amplitude formula is a general formula when Amp becomes 
15 higher in the low dot density area and lower in the high dot 
density area in Q value quantize processing, and the amplitude 
can be easily controlled by one formula. 

The meaning of this formula will be described with reference 
to Fig. 7 to Fig. 10. Fig. 7A shows the relationship between 
20 the average value Avg and the remainder mod when binary processing 
(that is, Q = 2) is performed. Fig. 8A shows the relationship 
between the average value Avg and mod (Avg, 255/ (Q-l ) )-255/2 
(Q-l). When a power method is used, the formula indicates the 
relationship between the average value Avg and the power value 
25 (mod (Avg, 255/ (Q-l)) - 255/2 (Q-l)) y, as shown in Fig. 9A. 
Therefore, the amplitude which is to be the high amplitude in 
the low dot density area and the amplitude which is to be the 
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low amplitude in the high dot density area are obtained from 
the average value (hereafter Avg) . In the formula ( 1 ) , the power 
value is divided by (S 4 + l ) , which is in proportion to the difference 
6. As Fig. 10 shows, the amplitude Amp in the formula (1) is 
5 high in a part where the density variation is low, since 5 is 
small. In other words, the function of the sine wave is robust, 
and the result of dot dispersion is strongly expressed. In a 
part where the density variation is high, the difference 6 
increases, and the amplitude Amp in the formula (1) decreases. 
10 Therefore, by adding the AM-FM screen function, unclear edges 

J;* can be prevented, and dots can be arranged more uniformly. 

;;S Fig. 7B shows the relationship between the average value 

•Li i 

^ Avg and the remainder mod when quaternary processing (that is 

'" : Q = 4 ) is performed. Fig. 8B shows the relationship between the 

;!* 15 average value Avg and mod (Avg, 255/ (Q-l)) - 255/2 (Q-l). When 

jj: a power method is used, the formula indicates the relationship 

]:=; between the average value Avg and the power value (mod (Avg, 

255/ (Q-l)) - 255/2 (Q-l)) y as shown in Fig. 9B. Therefore, the 
amplitude which is to be the high amplitude in the low dot density 

20 area and the amplitude which is to be the low amplitude in the 
high dot density area are obtained from the average value 
(hereafter Avg) in the same way. In the formula (1), the power 
value is divided by ( 5 4 + 1 ) , which is in proportion to the difference 
5. As Fig. 10 shows, the amplitude Amp in the formula (1) is 

25 high in a part when density variation is low, since 6 is small. 
In other words, the function of the sine wave is robust, and 
the result of dot dispersion is strongly expressed. In a part 
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where density variation is high, the difference 6 increases , and 
the amplitude Amp in the formula (1) decreases. 

Now the frequency Frq of the AM-FM screen function will 
be described with reference to Fig. 11 to Fig. 14. The frequency 
5 also depends on the input pixel value obtained by the detection 
circuit 7. The frequency is strongly correlated with the ideal 
dot interval in the target pixel value. For example, as Fig. 
11 shows, the input value of 16 x 16 pixels with uniform density 
level 247 is binarized by error diffusion. In this case, there 

10 are 8 black dots in 256 pixels as a result of error diffusion. 

In one-dimensional terms, 8 in 256 pixels means that one 
black dot appears at every 32 pixels . In two-dimensions , a square 
root is found, that is, black dots are arranged at an interval 
of average [square root sign] 32 = 4 [square root sign] 2 pixels, 

15 and this is an ideal interval of black dots in a two-dimensional 
space. 

Considering that a 45° dot arrangement is most comfortable 
for humans to view, dots are arranged in this way, then the dot 
interval in the horizontal and the vertical direction is 4 [ square 
20 root sign] 2 x [square root sign] 1/2x2 = 8, and one black 
dot is arranged at every 8 pixels. 

In other words, in the level 24 7, the cycle (1/Frg) of the 

AM-FM screen function is 2ji/8 = jc/4 , and at quantization processing 
black dots are trapped according to this cycle. Since black dots 
25 are dispersed on a white background in the highlight part, the 
control of black dots should be a primary consideration, and 
since white dots are dispersed on a black background in the shadow 
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part, the control of white dots should be a primary consideration 
as well. 

Given these considerations , frequency is used depending 
on the input pixel value as in the following general formula 
in a C language program. 

If mod (Avg, 255/ ( Q-l )) -255/2 ( Q-l )>=0 

Frq = [square root sign] (1/2) * [square root 

sign] [ [255/ (Q-l) ] / [ 255/ (Q-l ) -mod( Avg, 255/ (Q-l) ) ] ]/7T 

(2) 

else 

Frq = [square root sign] (1/2) • [square root sign] 
[ [ (255/(Q-l) ]/mod (Avg, 255/(Q-l))]/ K (3) 

Fig. 12 shows the relationship between the average value 
and the frequency in the case of binary processing, and in the 
case of quaternary processing. In other words, high frequency 
is set in the high dot density area, and low frequency is set 
in the low dot density area. As a result, the cycle of the AM-FM 
screen function becomes 1/Frg. 

Here the cycle of the AM-FM screen function is considered 
two-dimensionally . If the cycle in the horizontal direction is 
cx (= 1/Frq x ), and the cycle in the vertical direction is cy ( = 
1/Frq y ), then the AM-FM screen function (AF (x, y)) is given by 
the following formula. 

AF(x, y) =Amp • sin(cx + cy + 5) (6 : initial phase) (4) 
Therefore the AM-FM screen function indicates the high 
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amplitude and the low frequency in the low dot density area, 
and the low amplitude and the high frequency in the high dot 
density area. If only cx is incremented in one pixel unit at 
the same time as error diffusion processing, thenAF (x, y) becomes 
5 a one-dimensional oscillation, as shown in Fig. 13, and changes 
in the main scanning direction. 

In this state, however, a component of the sine wave in 
the vertical direction disperses. So two-dimensional control 
is required. For this, cy must also be incremented. However, 
10 if cx and cy are incremented simultaneously, AF (x, y) will be 
~ different from the intended waveform. In order to appropriately 

t control the AM-FM screen function in a two-dimensional space, 

7 the horizontal increment and the vertical increment are 

alternately switched so as to perform pseudo-two-dimensional 
* 15 control. 

II In other words, the switching conditions are to switch at 

i a predetermined pixel interval, and to switch when the edges 

of the image in the vertical direction are generated, as shown 
in Fig. 13. When these conditions are met, the progress of the 
20 sine wave is reset , cy is incremented, and the sine wave progresses 
in the vertical direction . To perform this switching process ing , 
a cycle of a sine wave of one previous scanning line recorded 
in the memory circuit 9 in Fig. 2 is used. 

As a result of these two switching conditions, the AM-FM 
25 screen function behaves two-dimensionally according to the 
density variation of the image. As Fig. 14 shows, the AM-FM screen 
function becomes the high amplitude and the low frequency in 
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the low dot density area, and the low amplitude and the high 
frequency in the high dot density area. 

Fig. 15 is a flowchart depicting half tone processing 
according to an embodiment of the present invention to perform 
5 the control. 

(51) Addresses x and y of the target pixel are initialized 
to " 1 " . The cycles cx and cy of the screen function are initialized 
to " 0 " . 

(52) Whether y is the number of vertical pixels or less 
10 is judged. Processing ends if y exceeds the number of vertical 

pixels . 

(53) Whether x is the number of horizontal pixels or less 
is judged. Processing advances to step SI 1 if x exceeds the number 
of horizontal pixels. 

15 (S4) Amplitude Amp is calculated from the above formula 

(1), and the frequency freq is calculated from the formulas (2) 
and ( 3 ) . 

(55) To judge the above mentioned switching conditions in 
the vertical direction, mod (x, cl) is calculated, and whether 

20 the calculated value is "0" is judged, cl is the above mentioned 
predetermined interval in the horizontal direction. Also, 
whether the difference 6 is 0 is judged. 

(56) If neither is "0", switching is not performed, so that 

jt/freq is added to the cycle cx in the formula (4). 
25 (S7) If either is " 0" , the above mentioned switching 

conditions are met, so switching to the vertical direction is 
performed. In other words, cx is set to "0", progress in the 
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horizontal direction is stopped, and the cycle rireki (x) of 
the sine wave of one previous scanning is set to cy. 

(58) Then the AM-FM screen function value is calculated 
using the formula (4), and the calculated value along with the 

5 quantization errors E from the peripheral pixels are added to 
the target pixel value I. Then quantization processing is 
performed and the output 0 is calculated. 

(59) The quantization errors E, generated at that time, 

are calculated and recorded in the error buffer. Also jt/freq 
10 is stored in the cycle rireki (x). 

(510) x is incremented and processing returns to step S3. 

(511) y is incremented and [processing] returns to step 
S2. By the error diffusion and the AM-FM screen function in the 
above serial processing, the dot pattern of the output image 

15 can be appropriately controlled. 

In the case of a color image, the initial phases (6 in the 
formula (4)) of the sine wave in each C, M and Y plane are set 

to different values, 61, 62 and 63, as shown in Fig. 16, so that 
the primary color dots of C, M and Y can be arranged together 

20 at relatively equal intervals. As Fig. 17 shows, arranging C, 
M and Y dots together increases the spatial frequency components 
of the half tone image, graininess decreases as a result, and 
the color gamut expands by mixing the colors of dots. 

Embodiment will be explained next. Fig. 18 is a diagram 

25 depicting an output pattern of a comparison example using the 
Floyd error diffusion method, Fig. 19 is a diagram depicting 
an output pattern of a comparison example when two-way scan and 
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error filter switching are used , and Fig. 2 0 is a diagram depicting 
an output pattern when the AM-FM screen function is added to 
the case in Fig. 19 according to the present invention. In the 
comparison example in Fig. 19 , worm at the left and fingerprint 
5 at the right in Fig. 18 are improved compared with the comparison 
example in Fig. 18 . In the output image according to the present 
invention in Fig. 20, dot patterns with less graininess, quite 
similar to half tone dots, can be obtained in the low dot density 
area compared with the dot patterns obtained by the comparison 

10 example in Fig. 19. 

The present invention has been described with the embodiment 
above, but various variant forms are possible within the scope 
of the present invention, and these are not excluded from the 
scope of the present invention. 

15 In the present invention, the amplitude and the frequency 

of the s ine wave to be superimposed are modulated us ing the average 
value of the target pixel value and the peripheral pixel values. 
Since the average value indicates high/low dot density, a sine 
wave with the amplitude and the frequency according to the dot 

20 density can be superimposed. Therefore, even if errors are 
diffused to maintain grayscale levels, a dot pattern which is 
uniform and cyclic, such as half tone dots, can be arranged in 
a low dot density area, and a high frequency dot pattern, which 
is closer to checkers, can be arranged in a high dot density 

25 area, and as a result, a visually fine half tone image can be 
obtained. 
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